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Abstract 

We rewrite the Zwanziger horizon condition in terms of the Kugo-Ojima 
parameter for color confinement. This enables one to explain which value of the 
Kugo-Ojima parameter is allowed if the horizon condition is imposed. Although 
all the calculations are performed in the limit of vanishing Gribov parameter 
for simplicity, the obtained value is consistent with the result of numerical 
simulations. Consequently, the ghost propagator behaves like free and the 
gluon propagator is non-vanishing at low momenta, in harmony with recent 
lattice results and decoupling solution of the Schwinger-Dyson equation. The 
Kugo-Ojima criterion is realized only when the restriction is removed. 
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1 Introduction 



To solve color confinement problem is equal to answering a fundamental question of 
how the non-Abelian gauge theory can be quantized in a non-perturbative manner 
beyond the perturbation theory, as represented by the Gribov problem [1] . This topic 
is currently under extensive studies due to a fact that color confinement is closely 
related to the asymptotic behavior of the ghost and gluon propagators in the deep 
infrared region, causing a challenge for numerical simulations on larger lattices. 

It was shown [2] that the partition function of the D-dimensional Euclidean Yang- 
Mills theory restricted to the first Gribov region for avoiding Gribov copies can be 
written in the form: 

:= y r'^5((9^j4)detMexp{-^yM + 7 / (fxh{x)\, (1.1) 

where Syu is the Yang-Mills action, M is the Faddeev- Popov operator M := —d^D^ = 
~9fj,{d^ + g-^pi^) and h{x) = h[£/]{x) is the Zwanziger horizon function given by 

h{x) := -/ d^ygf^^^^^^{x){M-Y''{x,y)9f^^''K^{y)- (1-2) 

Here the parameter 7 called the Gribov parameter is determined by solving a gap 
equation, commonly called the horizon condition: 

{h{x)y = {N^ - 1)D. (1.3) 

The action corresponding to the partition function (11. ip contains the non- /oca/ horizon 
term: 

/ d^xhix) := -/ d^x I d^ygf^^^^fix)iM-'f^ix,y)gf^^^^^^iy). (1.4) 

Later, it has been shown [Hill] that the non-local action (II. ip can be put in an 
equivalent local form by introducing a set of complex conjugate commuting variables 
and anticommuting ones, which is called the Gribov-Zwanziger (GZ) action. We 
do not write it explicitly, since we do not use it in this paper. The GZ theory is 
renormalizable to all orders of perturbation theory. Hence, the restriction to the 
(first) Gribov region Q makes perfect sense at the quantum level, and finite results 
are obtained consistent with the renormalization group. 

Note that Q is obtained as the local minima of the gauge-fixing functional 

F:=J d^x<(xX(x). (1.5) 

However, it is known that Q still contain Gribov copies, since there exist many local 
minima starting from the same The gauge field configurations that are absolute 
minima of the gauge-fixing functional F are known as the fundamental modular region 
(FMR) KpMR- However, it is extremely difficult to deal with the Yang-Mills theory 
by discriminating between A^mr and Q. 

The Kugo-Ojima color confinement criterion [51E1E], i-e., a sufficient condition for 
color confinement, is given by 

M^^(O) = 5^^m(0) = -5^^, t.e.,u{0) = -1, (1.6) 
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where the so-called Kugo-Ojima parameter u{0) is defined by the infrared (IR) limit 
/c^ ^ of the two-point function of composite operators: 

- ^] S^'^ik') = I d^a;e^^(--^)(0|T[(D,^)^(x)(^?^. x ^1^(2/)]|0). (1.7) 



However, it should be noted that the Kugo-Ojima color confinement criterion was 
obtained in the framework of the BRST quantization for the usual Faddeev-Popov 
approach, which corresponds to the 7 = case of the above Gribov-Zwanziger for- 
mulation. In this paper, we consider the D-dimensional Yang-Mills theory in the 
covariant gauge with a gauge fixing parameter a defined by 



where 



Z= J[d£/][d^][d'^][d'^]exp{iS'^ij}, (1. 



Qtot c I C 

^YM - — ^YM + ^GF+FP, 

Sym :=- I d^'x^^^, ■ 



Sgf+fp := I d^'xl^- + ■ ^ + 2^ ■ d^D^^ 



a 
2' 



D^'rf:={d^ + gs^^x)^, (1.9) 

and the dot and the cross are defined as 

^.^:=^^^^, (^ X i^)^ := /^^^^^^^, (1.10) 

using the structure constant of the gauge group G = SU{N). The Landau gauge 
corresponds to a = 0. We assume that the Euclidean result is obtained by the Wick 
rotation of the Minkowski one. 

The usual Faddeev-Popov approach does not take care of the Gribov copy problem. 
Therefore, if one begins to avoid the Gribov copy by restricting the space of gauge field 
configuration, it may happen that the Kugo-Ojima criterion u{0) = —1 based on the 
Faddeev-Popov approach does not necessarily hold. In fact, the direct measurements 
on a lattice [8] confirm u{0) = —0.6 ~ —0.8. To my knowledge, there is no theoretical 
explanation for this result. 

In this paper we discuss a relationship between the Zwanziger horizon condition 
and the Kugo-Ojima color confinement criterion. In other words, we study how 
the Zwanziger horizon condition imposes a constraint on the possible value of the 
Kugo-Ojima parameter. For this purpose, we give an estimation of the average of 
the horizon function {h{x))'^ in terms of the Kugo-Ojima parameter u{0). We obtain 
some relations connecting three body and four body two-point functions for gluon and 
ghost fields, which able to give somewhat stronger results than those obtained from 
the Schwinger-Dyson equations for the ghost propagator and the ghost-antighost- 
gluon vertex function. 

However, we restrict the calculation to the 7 = limit, i.e., in the usual Faddeev- 
Popov approach on which the Kugo-Ojima is actually based. The 7 7^ case will be 
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treated in a subsequent paper. Instead, various relations are obtained at arbitrary 
momentum k and the deep infrared (IR) limit k will be taken only in the final 
stage. Finally, we give an argument supporting that the result obtained in the limit 
7 = is not far from the correct value. 

2 Two-point functions of composite operators in a 
covariant gauge 

In this paper, we consider the two-point function for composite operators. We use an 
abbreviated notation for the Fourier transform: 

{<Pt<p2)k ■■= I rfW=(--^)(0|T[</)^(x)</.f (y)]|0). (2.1) 

In particular, we pay attention to the deep infrared limit A; — > 0. We begin with 
reproducing the result of Kugo |B] • 

First, we consider the three-body two-point function, {{gs^^ x Using 
the definition of the covariant derivative, 

{D^^)\x) := d,^\x) + {g^, x ^)\x) = d,^\x) + ^7/^^^i<(a;)<^^(a;), (2.2) 

we obtain 

((^< X <^) =((D^^)^^^)fc - 

(^^^'^ + (^^^"^);.) , (2.3) 
where we have used in the second equality the integration by parts 

dl ^ -ik^, (2.4) 

and a well-known relation [^lE] 

((Z}^^)^^^)fc = if^b^^. (2.5) 

It is very useful to introduce the one-particle irreducible (IPI) part by 

((^7=< X ^)^^^)f ^ := ((^7=< X ^)^^^)fc/(^V^)fc. (2.6) 

Similarly, we can define another IPI part: 

{^\gs^. X ^ = {^\gs^, X ^)^),/(^^^"^),. (2.7) 

Then the above result (12. 3p is also written in the form 

((^7< X ^)V"^)r = -^K (-1 + ^(^^^"'')."') • (2.8) 
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This result f l2.3p or fl2.8p is stronger than the result of the naive analysis of the 
SD equation for the ghost propagator, since only the contracted form k^{{g£^^ x 
appears in the SD equation for the ghost propagator in momentum space: 



= - X (2.9a) 

or 

^ = -fc^^^^ - ik^{{g£^^ X ^)^^^)^^^. (2.9b) 

In fact, fl2.3p or (12.80 automatically satisfies this SD equation. 

Next, we consider the four-body two-point function, {{g^^fj. x '^)'^[g.G/y x '^)^)k- 
It is also well known [HIE] that the following identify holds. 

= {{d^D,^)\gs^, X <^)^)fc, (2.10) 

or 

= {{8^-8^^)%^ X + {8^{g£^, X ^)\g^, x ^)^),. (2.11) 

Therefore, we have 

zF((^i4 X ^)^(^^/, X <^)^)fe = -A;2(<^^(^< X ^)^)fc. (2.12) 
We now define the modified IPI (mlPI) part of ((^^^ x ^)^(^^^ x ^)^)fc by 
((^7=< X ^)\g^, X 

: = ((^7^^ X '^)\gK X ^)^), - {{g^, x ^)^^^)r(^''^'')fcC^''(^7=^<^ x <^)^)r- 

(2.13) 

The IPI part should be defined from the connected part. Hence, in the above defi- 
nition, {{gs^^ X '^)^{gs^u X '^)^)fc must be replaced by ((^^^ x ^)^(^^/^ x <^)^)^°"'^ := 

{{g^, X <^)^((7< X - ^ ^gfABCgfADE^^B^D^^_^^^C^E^^_ That is tO 

say. 



((^< X ^)^(<7< X ^i^)r" 

={{g^, X '^)^^"^)rc^'''^"'')fc(^''(^< X + {{gK X ^)^(^^. X ^f)r. 

(2.14) 

Therefore, IPI and mlPI is related as 

{{g^, X ^)^(^7< X 

(2.15) 

In what follows, we use ml PI part rather than IPI one. 
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By using fl2.12p and (12.31) . then, it turns out that 

= - X 'rff)k - ik^iigK X X ^ 

= C<r^(^^. X ^f)r, (2.16) 

where we have used the definition of ((fif^ x 'lo)^^^)^^ and (^^(gfia^ x ■ 
Thus we have obtained an important relationship connecting the three-body two-point 
function and the four-body one: 

{^\gs^. X = ik,{{gs^, X ^)\g^, x ■ (2-17) 

This relation was already pointed out to hold in Kugo [B] by a diagrammatical con- 
sideration (without the explicit derivation). This relation should be compared with 

(EH. 

In order to extract more property of the IPI part of the four-body function, we 
note that it is also written as 

{{g.i^, X ^)^{gK X 
={{D^^)\gK X - {{D^^)^^'^),{^''{gK x ^^Y^ . (2-18) 

since the definition f l2.13p leads to 

{{gs^, X ^)^((7< X 
= {{D^^)\gK X ^)^)fc - {{d,^)\g< X ^)^)fc 

- {{gK X '^)^'^'')k{'^''{gK X ^f)\P' 

= {{D,^)^{gK X ^)^)fc - {{d,^)^{gK X <^)^)fc 

- ((D^<^) V^),(^^(^7^, X <^)^)f ^ + ((9^^)^^^)fc(^^(^?< X 
= {{D,^)^{gs^, X ^)^)fc - {-ik^){^^{gK x ^)^)fe 

- ((Z}^<^) V^),(^^(^7< X <^)^)i^^ + (-zA;^)(^^^^),C^^(^7i/. X ^)^)r- 

(2.19) 

By using (12.50 and (12.170 . therefore, (I2.18P is cast into 

{{D^^)\gs^, X ^)^),. 
= {{9^, X ^)^(^?< X ^")^)-i^^ + {{D,^)^^'')k{^^{gK X 

=5/((^7^P X ^)\gK X ^1^)^''' + ^Y^^kP{{9K x ^)^(<7< x '^1^)^''' 
= - ^) ((^-^P X X (2.20) 

Thanks to (12.100 . we can introduce the function u{kP') according to Kugo & Ojima [5]: 
{{D,^)\gK X ^r), := - ^) 5^^n(A;2), (2.21) 



or combining this with (12.51) yields 

{{D,<^nD^^f), = l^g,^ - 5^Mk') - ^S^^. (2.22) 

If we combine ( ^^M with flMTl) . therefore, {{g^/^ x 'rfyig^/^ x 'rf)^)'^^^^ must have 
the form: 

{{g^, X <^)^(^7< X = [^?^.n(A;2) + A;,A;,t;(A;2)]5^^, (2.23) 
where w(A;^) is an arbitrary function and a different parameterization 

w{k) := A;2^;(A;2), (2.24) 

reads 



{{g^, X ^)^(^7< X 



B\mlPI 



(2.25) 



3 Preparing the horizon function 



The horizon function is defined only in the Landau gauge. However, we collect the 
necessary information for the same type of function as the horizon function for arbi- 
trary covariant gauge. 

The vacuum expectation value (VEV) of the would-be horizon function is 



{h{x)) = -J rf^y(^?r^^<(x)(M-i)^^(x,y)^7r^X''(^)) (3-la) 

= - / d^y(^7r^^<(x)<^^(x)^?r^X''(2/)^"''(^)) (3-lb) 

= -/ d^'yiig^, X ^)^(x)((7^^ x '^)^(y)), (3.1c) 

where we have used the identity which holds for any functional f{A) of A: 

{f{A)^\x)^^{y)) = {fiA){ArY^{x,y)). (3.2) 
Therefore, the average of the VEV of the would-be horizon function reads 

V^' f d^x{h{x)) =Vo' f d^x [ d^'viig^, X ^)^{x){g^, x ^)^(y)) 



X ^)''{gs^, X ^)^)fc=o, (3.3) 

where Vd is the volume of the D-dimensional spacetime. 
The use of (ElQ]) and fl2Tfl) yields 

((^< X '^)^{g^, X 
= {{D,^)\gs^, X ^^)^)fc - {{d,^)\gs^, x '#)^), 

= - ^ - k,k''{^^)^ {{gs^, X ^)^(^7< X (3.4) 

= f^^. _ _ kM"^"^)^ n(A;2)5^^ - A;^A;,(^'^),A;2|;(A;2)<5^^, (3.5) 
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where we have used fl2.23p in the last equahty and we assumed that color symmetry 
is not broken so that 

{^^^^), = (3.6) 
Performing the Lorentz and color contractions, we obtain 

(((?i4 X ^Y{gs^^ X ^)^), 
= (A^2 - \)\(p - 1 - A;2(<^^)fe) u{k^) - A;2(^<^)feA;\;(A;2)] 
= (Ar2 - 1) \{p - \)u{}^) - k^l,'€'i)u{u{k^) + k^ie)"^. (3.7) 

Therefore we arrive at the average of the would-be horizon function written in terms 
of the Kugo-Ojima parameter and the ghost propagator at A; = 0: 



lim(((7=<x^)^((?i4x'^)^)fc 



= -{N' -\)\{D - l)n(O) - lim [r [n(0) + w{^)\ , (3.8) 

L fe2— >o J 

where we have used limfc2_^o[^^'^'(^^)] = limfc2_^o['W^(^^)] = w(0). 

It is observed that the dressing function Z{k'^') := —k'^{'io^)k of the ghost, 

fc2^o^ ^ m(0)+w(0) ' ^ ' 

is finite, namely, limfc2^o[^^(^^)/t]~^ 7^ 0, as long as {h(0)) is finite, provided 
that m(0) + w(0) ^ 0. 

Note that w{0) ^ means the existence of the massless pole in the second part 
of ((fif^^ X '£)^{g£^y X "iY)"^^^^ in fl2.25p . Usually, the absence of the massless pole 
in assumed w(0) = 0. However, this is not quite rigorous. Therefore, we leave a 
possibility w(0) 7^ in what follows. 

Up to here, all the equations hold for any covariant gauge, i.e., any value of the 
gauge fixing parameter a. In the Landau gauge a = 0, we can study the relationship 
between the horizon function and the ghost propagator in more detail, as shown 
below. 



4 Restricting to the Landau gauge 

In what follows, we impose the Landau condition a = 0. Then we have an additional 
identity: 

((^7< X <^) V^)r'°=° = ((^7< X ^)\gs^^ X ^)^)f ^^'-=0(-zA;.). (4.I) 

This is because the derivative factor acting on the antighost at the ghost-gluon- 
antighost vertex at the right end of the diagram can be transferred to act on the 
external ghost, since (9^^ = in the Landau gauge. 
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Substituting fl2.3p and fl2.23p into fl4.ip . we obtain a relationship between the 
Kugo-Ojima function and the ghost propagator in the Landau gauge: 

Z(fc2)"=o := -A;2(^V^)^=o = [1 + u{k^) + ev{e)]-^6^^. (4.2) 

In the Landau gauge, therefore, we find using (13. 7p 

(((?<, X ^^g^, X ^)^)r° 
= {N^ - 1) {(D - 1) u{e) - A;2(^^)-=0[u(A;2) + kMk% 

and the average of the horizon function reads 

(/i(0))"=o =1^^'/ rf^x(/i(x))"=° 

= - hm {{g^, X ^)^(^?^^ X 

= -iN'- 1) ((/} - 1)^(0) - \uay{^'^)r'MO) + w 

--(iv--i) (fl-iHo)+ (4.4) 

[ 1 + u(0) + w{0) J 

The same result is obtained even if we start from the original definition (I2.13P by 
taking into account (12.170 and (14.10 : 

{{g^, X ^)''{gK X 
■. = {i9K X ^)^(<7< X + (((7^^ X ^)^^^y/'{^^^^),{^''{gK X 

=((^^M X ^)^igK X 

+ k^k.iig^^ X ^)^((7< X ^^(^^^^)fc(((7^p X ^)''{gK x 

= [^?M^w(A;') + k,kMk')]S^'' + A;^A;,(^^^^)fc[w(A:2) + kMk'W, (4.5) 

where we have used (I2.23P in the last step. After performing Lorentz and color 
contraction, indeed, we obtain the same result as (14.31) : 

{{g^, X ^f{gs^, X ^)^)r' 
= {N'' - 1) {Dui^) + ev{k^) + e{^'i)k[u{k^) + k^v{k^)f] 

Note that the horizon function (14.41) defined above differs from that obtained by 
taking the IR limit of ^TF^: 

- \^^{{D,^)\D,^)^)T' = -{N' - 1) {{D - IHO) - 1} , (4.7) 

except for a special value of m(0) (and w{0)). This is because the total derivatives 
in the Fourier transform can not be discarded in the above calculations. Indeed, we 



8 



have taken into account all the total derivative terms. Otherwise, i \4.7\i agreed with 

(531)0 

In the Landau gauge, we find 

= - z{{gs^^ X ^^y{gK X ^1^)^''''"=°^., (4.8) 
which is contracted with —ikf^ to give 

= - k,{ig^, X X ^f)f''''"=\,. (4.9) 

Then the inverse ghost propagator obeys 

(^^^^),-^ = - e6^^ - k.iig^, X ^)^igK X ^)^)r''^'"=°fc, 

= - fc^n;^,^(fc)fc., (4.10) 

where we have introduced 

K^ik) :=g,J^^ + {{g^, x ^)^(^< x ^-)^)™i^^'-o 

=5^^{[1 + u{e)]g,, + k.kMk')}- (4.11) 

Thus we have another relationship between the IR limit of the ghost dressing function 
and the Kugo-Ojima parameter (up to the value of w{0)): 

\un\-e{^^^%]-' = \im\k,U^^^^{k)k^/e] 

=5^^[l + n(0) + u; 



nr(o)/^ + (1-1/^)^(0). (4.12) 



5 Horizon condition and Kugo-Ojima parameter 

The horizon condition is written as 

{h{0)y = {N^ -1)D. (5.1) 

The actual value 7* of the Gribov parameter 7 is determined by solving this gap 
equation in a self-consistent way. Suppose that 7 is so small that the left-hand side 
can be expanded in powers of 7 around 7 = 0: 

Vd 97 



1 d\nZ^ 



+ 7 ^ 



7=0 



Vd 972 



+ 0(7^ 

7=0 



(/z(0)),=o + 7 / d^'yihiO); h{y))~'- + 0(f). (5.2) 
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Figure 1: The plot of {h{0))'^ ^ versus u{0) for = 4. 



Here the expectation value at 7 = is to be calculated in the usual Faddeev-Popov 
approach. In fact, we have obtained {h{0))^=Q as a function of n(0) in fl4.4p . 

We consider how the Zwanziger horizon condition (of reducing the integration 
region to the first Gribov region) restricts the possible value of the Kugo-Ojima pa- 
rameter m(0). If 7* is small, the horizon condition is approximately written by using 
(/i(0))?=0- 



'7= 



imns = -{N'- 1) {(D - i)u(o) + ^ (N^- - m (5.3) 

First, we consider the case of w(0) = 0. Then the horizon condition at 7 = is 

(M0))^:o° = - (iV^ - 1) \{D - 1)^(0) + Y^^} = - 1)^- (5-4) 

See Fig. □ for a plot of (/i(0))^=o versus u(0) for D = A. Note that n(0) 
in the vanishing coupling limit (7^0. Therefore, the allowed region for u(0) is 
u(0) e (—1, 0] for (/i(0))"rQ G [0, +00] in a branch going through the origin u(0) = 0. 
The horizon condition (15.31) leads to the second degree algebraic equation for u(0), 
{D — l)u^ + 2Du + D = 0, whose solutions are given by m = {—D ± ^/D)/{D — 1) 
e.g., for D = 4, 3u^ + 8m + 4 = has solutions m(0) = —2/3, —2, and we adopt 

u{0) = -| (5.5) 

because u{0) = —2/3 belongs to a branch going through the origin, while this is not 
the case for m(0) = —2. 



^This should be compared with the claim of Dudal et al. [16] 
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Even for 7 = 0, the Zwanziger horizon condition (/i(0))"^° = (A^^ — 1)D does 
not agree with the Kugo-Ojima color confinement criterion n(0) = —1. In this ap- 
proximation 7* ^ 0, the Kugo-Ojima criterion n(0) = —1 is reahzed only when 
(/i(0))-y=o = +00, which implies that there is no restriction for the gauge field space. 

Within this approximation 7* ~ 0, the horizon condition forces the ghost prop- 
agator to behave like free l/k"^ aX k = 0, no more singular than for D = A, 
irrespective of the number of color A^, 

lim [-A;2(^^^^)fc]-i = + u{0)] = ^5^^ ^ 0, Z{0) = 3. (5.6) 

Incidentally, the direct measurements of the Kugo-Ojima parameter u{0) on a 
lattice have been performed by imposing the absolute Landau gauge which restricts 
the integration region to the FMR. The result [8] is u{0) ^ —0.6 ~ —0.8. Our rough 
estimate u{0) = —2/3 = —0.67 is very near to this result of numerical simulations. 
This suggests that 7* is relatively small. See [13] for another support for this approx- 
imation. 

Incidentally, the formal power series, j^^^l^j = ^(0)[1 + u{0)]~^ = u{0) — u(0)^ + 
u{0)^ + ■ ■ ■ yields the horizon condition 

(/^(0))-r° =(iV2 - 1) {-DuiO) + n(0)2 - u{Of + ...]^{N^- l)D. (5.7) 

If we took into account only a linear term in m(0) on the left-hand side, then the 
Kugo-Ojima criterion n(0) = —1 would be satisfied and the ghost dressing function 
Z(0) = [1 + n(0)]-i would diverge. 

From the viewpoint of perturbation theory in the coupling constant (7, the Kugo- 
Ojima parameter m(0) begins with the order g'^, i.e., u(0) = uig"^ + U2g'^ + ■ ■ ■ . If we 
compare both sides of fl5.7p order by order in the coupling constant, then the Kugo- 
Ojima criterion is satisfied u(0) = —1 in the lowest 0{g^), and the ghost dressing 
function Z{0) = [l+u{0)]-^ diver ges, in agreement with the Gribov original result [1] 
where 0{g^) terms are taken into account. 

Second, we consider the case of w(0) 7^ 0. The horizon condition alone is not 
sufficient to determine both m(0) and w(0). Suppose the Kugo-Ojima confinement 
criterion is satisfied m(0) = —1. Then the horizon condition at 7 = is 

m)-Z^ = -{N'- 1) |-D + 1 + ^^^^^} = (N' - 1)D. (5.8) 

This leads to the value of w{0) irrespective of the spacetime dimension D and the 
number of color A^: 

w{0) = 1/2 for any D. (5.9) 

Even if m(0) = —1, therefore, the ghost propagator behaves like free 1/fc^ at /c = 0, no 
more singular than l/k"^: irrespective of the spacetime dimension D and the number 
of color 

lim [-A;2('^^<^^)fc]-i = d^^'wiO) = ^5^^ ^ 0, Z{0) = 2. (5.10) 
>o 2 

Thus, the ghost propagator behaves like free at low momenta, while the gluon 
propagator is non-vanishing at low momenta. 
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6 Conclusion and discussion 



In this paper, we have rewritten the Zwanziger horizon condition in terms of the 
Kugo-Ojima parameter u{0) for color confinement. Then we have examined the region 
allowed for the Kugo-Ojima parameter under the horizon condition. Although all the 
calculations are performed in the limit of vanishing Gribov parameter for simplicity, 
the obtained value u{0) = —2/3 for D = 4 is consistent with the result of numerical 
simulations [S]. 

Consequently, the ghost propagator behaves like free and the gluon propagator is 
non- vanishing at low momenta, in harmony with recent lattice results 0110], decou- 
pling solution of the Schwinger-Dyson equation [TT] and other approaches [T2 l[T3l[T^ . 

Our result suggests that the Kugo-Ojima criterion for color confinement u{0) = — 1 
will be realized only when the Zwanziger horizon condition is removed (or weakened) 
in a certain way. In order to clarify this issue, the 7 7^ case [15] will be treated in 
a subsequent paper. 

[Note added] In preparing this paper, an interesting paper by Zwanziger was 
posted to the archive [T7], in which the 7 7^ case was treated. The essential 
difference between the Zwanziger result for 7 7^ [T7] and ours for 7 = comes from 
the 2nd term ((c/j^ x '^)^^^)i^^('^^^^)fc(^^(c/^^ x which is written as 

^sing{k) in [l7]. At least in the 7 = case, namely, the usual FP approach, Asing{k) 
does not vanish and remains non- zero even after taking the k = limit, as we have 
examined in the above analyses. In fact, if we neglected the contribution from 

As^ngik) := X '^)^'^'')r('^'''^'')fcC^''(^7^. X ^f)r, (6.1) 

the average of the horizon function became instead of (15. 3p equal to 

(/i(0))"=o = ~{N^ - 1)Dm(0), (6.2) 

and the horizon condition (/i(0))°^° = (A^^ — 1)D would become equivalent to the 
Kugo-Ojima criterion n(0) = — 1. There may exist some discrepancy between 7 = 
and 7 7^ 0, if both papers are correct. 
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